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1. Introduction 



In an earlier paper |Jl|] we presented low-temperature series for the spontaneous 
magnetisation, susceptibility and specific heat of the spin-1 Ising model on the square 
lattice. In this paper we extend this work to higher spin values {S = |, 2, |, and 3). 
From general theoretical considerations, in particular renormalization group theory, it is 
expected that the critical exponents (at the physical singularity) depend only upon the 
dimensionality of the lattice and on the symmetry of the ordered state, and thus does 
not vary with spin magnitude. Numerical work on the Ising model with S* > 1 is quite 
sparse and little has been published since the mid-70's. Low-temperature expansions 
were obtained by Fox and Guttmann for 5* = 1 and S* = | for various two and three 
dimensional lattices. High-temperature expansions have been reported by a number 
of authors p, ^, |^ who mainly focused on three dimensional lattices. Generally the 
numerical work has confirmed spin independence. Recently, Matveev and Shrock 
studied the distribution of zeros of the partition function of the square lattice Ising 
model for 5 = 1, |, and 2. While the physical critical behaviour of the spin-S* Ising 
model is fairly well understood, little is known about the non-physical singularities. One 
major reason for seeking more knowledge about the complex-temperature behaviour is 
the hope that this will help in the search for exact expressions for thermodynamic 
quantities which have not yet been calculated exactly. 



2. Low temperature series expansions 

The Hamiltonian defining the spin-S Ising model in a homogeneous magnetic field h 
may be written: 

n = j-,Us''^^^,) + ^Us-^.) (1) 

where the spin variables cxj may take the (25" + 1) values ai = S, S — 1, . . . , —S. The 
first sum runs over all nearest neighbour pairs and the second sum over all sites. The 
constants are chosen so the ground state (cTj = S Vi) has zero energy. The low- 
temperature expansion, as described by Sykes and Gaunt |^, is based on perturbations 
from the ground state. The expansion is expressed in terms of the low-temperature 
variable u = exp(— /JJ/S"^) and the field variable /i = exp{—ph/S), where P = 1/kT. 
The expansion of the partition function in powers of u may be expressed as 



oo 
k=0 



(2) 
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where ^a:(/^) are polynomials in fi. It is more convenient to express the field dependence 
in terms of the variable x = 1 — ^ 



k=0 

Using the standard definitions, we find the spontaneous magnetisation 

1 dlnZ 



(3) 



M{u) = M(0) + 



p dh 



S + Z^{u)/Zoiu), 



(4) 



h=0 



since x = in zero field. For the zero-field susceptibility we find 



dM 



dh 



h=0 



d_ U_,dz^ 

dh \ dh 



h=0 



Z2{u) Ziiu) ( Ziiu) 



Zq{u) Zo{u) \Zo{u)^ 



(5) 



The specific heat series is derived from the zero field partition function (via the 
internal energy U = —{d/dp) InZo), 



dU d"^ 



2lnZo = (/5J/^2)MM— InZo(w). 



(6) 



Thus in order to calculate the specific heat, spontaneous magnetisation and 
susceptibility one need only calculate the first three moments (with respect to x), Z^iu) 
for /c < 2, of the partition function. These moments are most efficiently evaluated using 
the finite lattice method. The algorithm was described in an earlier paper [|I|. For our 
present purpose it suffices to note that the infinite lattice partition function Z can be 
approximated by a product of partition functions Z^^n on finite {m x n) lattices. 



HZ 



mn 



with m <n and m + n < 



(7) 



m,n 



The weights a^n were derived by Enting p|, and are modified in the present 
algorithm to utilize the rotational symmetry of the square lattice. The number of terms 
derived correctly with the finite lattice method is given by the power of the lowest-order 
connected graph not contained in any of the rectangles considered. We use the time- 
limited version of the algorithm in which the largest rectangles are determined by 
a cut-off parameter 6max, m + n < r = 36max + 2. The simplest connected graphs not 
contained in such rectangles are chains of r sites all in the '5* — 1' state. From (|l|) we see 
that such chains give rise to terms of order 2{r + l)[S'^ — S{S — l)] + {r — l)[S'^ — {S — lY] = 
r{AS — 1) + 1, from the 2(r + 1) interactions between spins in states 'S" and 'S* — 1' 
and the r — 1 interactions between spins both in state '5 — 1'. For a given value of 6max 
the series expansion is thus correct to order ■u(3Vax+2)(45-i)_ earlier paper []T| we 
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reported on the 5 = 1 case where we went to 6max = 8 giving a series correct to u'^^. 
We have since extended these series to u^^^ using a more efficient parallel algorithm and 
a new extrapolation procedure |p. For the present work we have calculated the series 
expansions for S' = |, 2, | and 3, deriving series correct to (&max = 6) for S" = |, 
u"' (&max = 5) for S = 2, u^^^ = 4) for 5 = |, and n^^^ {b^^^ = 4) for 5 = 3. 

3. Analysis of the series 

The series for the spontaneous magnetisation, the susceptibility and the specific heat of 
the spin-S* Ising model are expected to exhibit critical behaviour of the forms 



M{u) ~ — uy^^[l + cij^i{uj — u) + aj^/^{uj 



X{u) ~ n B,{u, - u) [1 + 6,- -u) + b,,^{u, - n)^^ + ...], (9) 



Cv{u) ~ W Cj{Uj - U) [1 + Cj,i{Uj -U)+ Cj^A{Uj - u)^' +...], (10) 

j 

where the terms involving represent the leading non-analytic confluent singularity 
and the dots . . . represent higher order analytic and non-analytic confluent terms. By 
universality it is expected that the leading critical exponents at the physical singularity, 
Uc, equal those of the spin-^ Ising model, i.e., = |, 7' = | and a' = (logarithmic 
divergence). 

We analysed the series using differential approximants (see Ref. |l^ for a 
comprehensive review), which allows one to locate the singularities and estimate the 
associated critical exponents fairly accurately, even in cases such as these where 
there are many singularities. We find that ordinary Dlog Fade approximants (ffist 
order homogeneous differential approximants) yield the most accurate estimates for 
the physical singularity of the magnetisation series, whereas ffist- and second-order 
inhomogeneous approximants are required in order to analyse the susceptibility and 
specific heat series. Here it suffices to say that a i^'th-order differential approximant 
to a function / is formed by matching the ffist series coefficients to an inhomogeneous 



differential equation of the form (see |jTO| for details) 



j2Q,{x){x^yf{x)=P{x) (11) 



where Qi and P are polynomials of order Ni and L, respectively. First and second order 
approximants are denoted by [L/Nq; Ni] and [L/Nq; Ni, N2], respectively. 
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3.1. The physical singularity 

In this section we focus on the behaviour at the physical critical point. First we give a 
somewhat detailed summary of the analysis of the spin-| series so as to introduce the 
various techniques and approximation procedures that we have applied in the analysis. 
Generally the estimates for the critical parameters at the physical singularity are quite 
poor because the series have many non-physical singularities closer to the origin and 
even for the spin-1 series |jl|, P] the convergence of the estimates to the true values of 
the critical parameters is very slow. We see no evidence that the critical exponents of 
spin-S* Ising model aren't in agreement with the universality hypothesis. Under this 
assumption, we have derived improved estimates for the location of the physical critical 
point and the critical amplitudes. 

In table |I] we have listed the estimates for the physical singularity and critical 
exponent for the spontaneous magnetisation of the spin-| Ising model. The estimates 
were obtained from homogeneous differential approximants (which are equivalent to 
Dlog Fade approximants). There is a quite substantial spread among the various 
approximants with most approximants yielding estimates around Uc — 0.7380 and 
P ~ 0.130. The estimates of (3, while generally on the large side, are consistent with 
expectations of universality which would indicate that /? = |. If we assume this value 
to be exact we see that the approximants (assuming a linear dependence of (3 on Uc) 
would lead to Uc — 0.73775. 

Table 1. Estimates for Uc and /3 for the spin-| Ising model as obtained from [A^, M] 
homogeneous first-order differential approximants. 



N 


[N-1,N] 


[N,N] 


[N+ 


1,N] 




Uc 


(3 


Uc 


/3 


Uc 




40 


0.738148 


0.1306 


0.738167 


0.1308 


0.738049 


0.1295 


41 


0.738124 


0.1303 


0.738020 


0.1291 


0.738081 


0.1298 


42 


0.737908 


0.1275 


0.737948 


0.1281 


0.737125 


0.1085 


43 


0.737918 


0.1277 


0.738046 


0.1294 


0.738099 


0.1300 


44 


0.738128 


0.1303 


0.738105 


0.1301 


0.738098 


0.1300 


45 


0.738123 


0.1303 


0.738059 


0.1296 


0.740267 


0.1038 


46 


0.737958 


0.1283 


0.738135 


0.1304 


0.738140 


0.1304 


47 


0.738140 


0.1304 


0.738135 


0.1304 


0.738331 


0.1317 


48 


0.736928 


0.1047 


0.737705 


0.1242 


0.737673 


0.1236 


49 


0.737676 


0.1236 


0.737700 


0.1241 


0.737867 


0.1271 


50 


0.738187 


0.1313 


0.737810 


0.1261 







In table |^ and table ^ we have listed estimates for the position of the physical 
singularities and critical exponents of the series for susceptibility and specific heat of 
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the spin- 1 model. Since the first non-zero term in these series is the estimates were 
obtained by analysing the series x('u)/m^ and Cy{u)/u^ . The estimates were obtained by 
averaging first order [L/N\ M] and second order [L/N; M; M] inhomogeneous differential 
approximants with \N — M\ < 1. For each order L of the inhomogeneous polynomial 
we averaged over most approximants to the series which as a minimum used all the 
series terms up to the last 15 or so. Some approximants were excluded from the 
averages because the estimates were obviously spurious. Examples include the [47,48] 
and [46,45] approximants in table 0. The error quoted for these estimates reflects the 
spread (basically one standard deviation) among the approximants. Note that these 
error bounds should not be viewed as a measure of the true error as they cannot include 
possible systematic sources of error. While the estimates aren't very good, we see 
that the estimates for Uc are consistent with the value Uc — 0.73775 obtained from the 
magnetisation series by demanding /9 = | and that the exponent estimates are consistent 
with universality expectations of 7' = 7/4 and a' = 0. 

Table 2. Estimates for Uc and 7' for the spin-| Ising model as obtained from 
inhomogeneous first- and second-order differential approximants. 





1. order DA 


2. order DA 


L 


Uc 


7' 


Uc 


7' 





0.73787(40) 


1.848(63) 


0.73802(37) 


1.864(58) 


1 


0.73808(3f) 


1.882(52) 


0.73810(26) 


1.868(49) 


2 


0.73800(19) 


1.864(34) 


0.73818(20) 


1.882(39) 


3 


0.73804(23) 


1.874(43) 


0.73804(33) 


1.848(72) 


4 


0.73792(48) 


1.82(10) 


0.73805(38) 


1.863(69) 


5 


0.73803(46) 


1.895(65) 


0.73808(25) 


1.852(69) 


6 


0.73787(53) 


1.839(99) 


0.73803(53) 


1.82(18) 


7 


0.73823(18) 


1.50(98) 


0.73792(51) 


1.80(12) 


8 


0.73774(64) 


1.76(20) 


0.73808(31) 


1.861(62) 



As for the critical exponents, it is obvious that the behaviour at Uc (except for S* = | 
and 1) isn't represented very well by the series. This discrepancy, which becomes more 
pronounced as 5* increases, is hardly surprising given that the number of non-physical 
singularities within the physical disc increases rapidly with spin magnitude (see following 
section for details). The quite complicated singularity structure of the series simply 
tends to obscure the behaviour at the physical singularity. This problem is possibly 
further aggravated by the presence of confluent terms. The only series which yields 
reasonably accurate estimates is the magnetisation from which we estimate /3 = 0.139(4), 
0.138(5) and 0.132(2) for 5 = 2, | and 3, respectively. Again, the quoted errors are 
merely a measure of the spread among the approximants rather than the true error. The 
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Table 3. Estimates for Uc and a' for the spin-| Ising model as obtained from 
inhomogeneous first- and second-order differential approximants. 



1. order DA 2. order DA 



L 




a' 


Uc 


a' 





0.74062(88) 


0.343(16) 


0.7393(18) 


0.17(25) 


f 


0.74030(88) 


0.320(80) 


0.7382(15) 


0.20(80) 


2 


0.7397(20) 


0.24(32) 


0.7389(18) 


0.12(20) 


3 


0.740f(f0) 


0.32(10) 


0.7384(16) 


0.07(23) 


4 


0.7370(28) 


0.06(71) 


0.7381(17) 


0.03(31) 


5 


0.738f(2f) 


0.04(38) 


0.7378(21) 


0.05(48) 


6 


0.7373(25) 


0.24(61) 


0.7388(29) 


0.07(53) 


7 


0.7357(24) 


0.21(64) 


0.7381(28) 


0.33(90) 


8 


0.7356(24) 


0.25(68) 


0.7386(25) 


0.02(66) 



differential approximant analysis of the higher S series for the susceptibility and specific 
heat yields little of value. Estimates for the critical exponent 7' fluctuate wildly and 
lie somewhere between 0.5 and 2 while generally favouring values below 7/4. Similarly, 
estimates for a' lie between -0.5 and 1. So while no sensible estimates can be obtained 
there is no evidence to suggest that the exponents are not consistent with universahty. 

While this situation is somewhat disappointing it is hardly surprising in light of the 
behaviour of the spin-1 series, where our earlier analysis showed a very slow convergence 
of estimates towards the true values of the critical parameters [|T], ^. Though the order 
to which the higher spin-S* series are correct exceeds that of the spin-1 series this is really 
just a consequence of the definition of the expansion variable u. We would expect the 
accuracy of estimates to depend not so much on the actual order of the series as much 
as on the maximal cut-off given by &max- In essence the accuracy is determined by the 
number of distinct graphs, consisting of spins flipped from the ground state (irrespective 
of the actual value of the spins), that one has summed over. One should therefore not 
expect more accurate estimates from the higher spin-S" series than those one could have 
obtained by truncating the spin-1 series at an order determined by the associated value 

of &max- 

One may hope to obtain improved estimates for Uc by raising the relevant series to 
the power 1/A, where A is the expected leading critical exponent, and look for simple 
zeros and poles of the resulting series. This procedure of biasing works quite well for the 
magnetisation and susceptibility series (it obviously cannot be used for the specific heat 
series). It is well known that the analysis of series exhibiting a logarithmic divergence, 
as we expect of the specific heat series, is particularly difficult. A fairly simple way of 
circumventing these problems is to study the derivative of the specific heat, d/duCv{u). 
The series for this quantity should have a simple pole at Uc, a situation much more 
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amenable to analysis by either differential approximants or even just ordinary Pade 
approximants. This approach does indeed confirm the logarithmic divergence at Uc, 
though the evidence becomes rather circumstantial for higher values of S. The estimates 
for Uc derived in this fashion are tabulated in table ^ and were obtained by averaging 
ordinary [A^ + K, N] Pade approximants {K = 0, ±1) with 2N + i^' + 15 not less than 
the order of the series. The error quoted for these estimates again merely reflects the 
spread among the approximants. 



Table 4. Biased estimates for the physical singularity. 
S Magnetisation Susceptibility Specific Heat 

f 0.73774(2) 0.7372(2) 0.7379(5) 

2 0.8293(2) 0.8288(2) 0.833(3) 
f 0.8795(3) 0.881(3) 0.882(2) 

3 0.9107(4) 0.914(1) 0.905(4) 



It is often possible to find a transformation of variable which will map the non- 
physical singularities outside the transformed physical disc. One such transformation is 
given hj u = x/{2 — x). Although the series in the transformed variable have radii of 
convergence determined by the physical singularity, this transformation turns out to be 
of little use and does not allow us obtain better estimates for the critical parameters. 
This is probably because there still are singularities close to the physical disc and because 



such singularity- moving transformations may introduce long-period oscillations |]T0 



We have calculated the critical amplitudes using two different methods, both of 
which are very simple and easy to implement. In the first method, we note that if 
f{u) ~ ^4(1 — u/uc)~^, then it follows that {uc — u)f^^^\ 

u=uc ~ A^^^Uc- So we simply 
form the series for g(u) = {uc — u)f^^^ and evaluate Pade approximants to this series 
at Uc- The result is just A^/^Uc- This procedure works well for the magnetisation and 
susceptibility series (it obviously cannot be used to analyse the specific heat series). For 
the specific heat series two different approaches have been used. In the first approach we 
use the 'trick' applied previously and look at the derivative of the specific heat series for 
which the above method should work with A = 1. In table § we have listed the estimates 
for the critical amplitudes obtained in this fashion. As usual estimates for any given 
value of Uc were obtained by averaging over many higher order approximants, and the 
error-estimates in table ^ reflect both the spread among the various approximants as well 
as the depence on Uc- In the second approach we start from f{u) ~ y41n(l — u/uc) and 
form the series g{u) = exp(— /(«)) which has a singularity at Uc with exponent A. One 
virtue of this approach is that no prior estimate of Uc is needed. However, the spread 
among estimates from different approximants is very substantial though consistent with 
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table ^. Biasing the estimates at Uc also confirms the value of the amplitude though 
generally the spread is larger than for the first approach. For the spin-3 susceptibility 
and specific heat series we could not obtain reliable amplitude estimates since the spread 
tended to be larger than the average value and the poor estimate of Uc leads to even 
greater errors. 

Table 5. Estimates for the amplitudes at the physical singularity. 



s 


Magnetisation 


Susceptibility 


Specific Heat 


3 
2 


1.875(5) 


0.019(3) 


52(2) 


2 


2.57(2) 


0.0088(5) 


110(5) 


5 
2 


3.33(3) 


0.006(2) 


190(10) 


3 


4.10(5) 







In the second method, proposed by Liu and Fisher |ri|, one starts from f{u) ~ 
A{u){l — u/uc)~'^ + B{u) and then forms the auxiliary function (/(u) = (1 — u/mc)^/(m) ~ 
A{u) + B{u){l — u/uc)^- Thus the required amplitude is now the background term 
in g{u), which can be obtained from inhomogeneous differential approximants p!0[ . 
This method can also be used to study the specific heat series. One now starts 
from f{u) ~ A{u) ln{l — u/uc) + B{u) and then looks at the auxiliary function 
g{u) = f{u)/ ln(l — u/uc). As before the amplitude can be obtained as the background 
term in g{u). This analysis yields amplitude estimates consistent with those in table |, 
but with larger error-bars. 

In table ^ we have listed our final estimates for the physical singularities and the 
associated exponents and amplitudes. For the estimates of the position of the physical 
singularities we have placed most weight on the biased analysis of the magnetisation 
series. In the spin-| case, Uc and the exponents a' and (3 and the amplitudes Ac and 
Am are known exactly due to the calculation of the free energy by Onsager [jl2| and 
the magnetisation by Yang |]T3| . The susceptibility amplitude is known to very high 
precision |]14|. The spin-1 estimates are from H. 



3.2. Non-physical singularities 

Except for S = \ the series have a radius of covergence smaller than Uc due to 
singularities in the complex w-plane closer to the origin than the physical critical point. 
Since all the coefficients in the expansion are real, complex singularities always come in 
pairs. The number of non-physical singularities appears to increase quite dramatically 
with S thus making it exceedingly hard to locate them accurately for large S. 

In order to locate the non-physical singularities in a systematic fashion we used the 
following procedure: We calculate all [L/N]M] inhomogeneous first order differential 
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Table 6. The physical singularities and associated exponents and amplitudes. 



s 


Uc 


P 


Am 


i 


Ax 


a' 


Ac 


1 

2 


3- 2^/2 


1 

8 


1.138789 


7 

4 


0.584850 





5.40658 


1 


0.5540653(5) 


0.12507(3) 


1.2083(2) 


1.750(1) 


0.0617(1) 


0.0005(10) 


22.3(5) 


3 
2 


0.73775(15) 


0.128(3) 


1.875(15) 


1.85(15) 


0.019(5) 


0.0(3) 


52(4) 


2 


0.8293(3) 


0.139(4) 


2.57(4) 




0.009(1) 




110(10) 


5 
2 


0.8795(5) 


0.138(5) 


3.33(6) 




0.006(2) 




190(20) 


3 


0.911(1) 


0.132(2) 


4.1(1) 











approximants with — M\ < 1 using all or almost all series terms for 10 < L < 16. 
(We discard no more than the last 15-20 terms.) Each approximant yields M possible 
singularities and associated exponents from the M zeroes of Qi (many of these are of 
course not actual singularities of the series but merely spurious zeros of Qi). Next we 
sort these 'singularities' into equivalence classes by the criterion that they lie at most a 
distance apart. An equivalence class is accepted as a singularity if it contains more 
than Na approximants {Na can be adjusted but we typically use a value around 2/3 of 
the total number of approximants), and an estimate for the singularity and exponent 
is obtained by averaging over the approximants (the spread among the approximants is 
also calculated). This calculation is then repeated for k — 1, k — 2, ... until a minimal 
value of 5 or so. To avoid outputting well converged singularities at every level, once 
an equivalence class has been accepted, the approximants which are members of it are 
removed, and the subsequent analysis is carried out on the remaining data only. This 
procedure is applied to each series in turn producing tables of possible singularities. 
Next we look at these tables in order to identify the true singularities. 

In table we have listed the non-physical singularities that we believe to have been 
identified with some degree of certainty and accuracy. For higher spin values several of 
these are marred by large error bounds and it is quite possible that we haven't been able 
to locate all non- physical singularities of the series, particularly for S' = | and 3. First 
we accepted any singularity which appeared in all the series at a reasonably early level, 
say k > 10. These singularities are marked 1 in table |^ and all of them are undoubtedly 
true singularities. Singularities which appear for A; < 10 are a lot more tricky to deal 
with. Generally we also expect that a singularity which appears for = 8 or 9 (or 
higher) in all series and for the majority of values of L is a true singularity of the series 
(these are marked 2 in table |^). However, we often find that some singularities appear 
for /c = 8 or higher in some series but at lower values of k all the way down to 5 in other 
series, and it is not easy to determine which ones are true singularities and which ones 
are not. Those marked 3 appear in all series and for all values of L while those marked 
4 appear in some series for all L but not neccesarily for all L in other series. 
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Table 7. Non-physical singularities Us and associated exponents of the spin-^* series. 





Us 


\Us\/Uc 


f3 


7' 


a' 


Spin-1 


1 


-0.598550(5) 


1.08 


0.1248(3) 


1.750(5) 


0.005(10) 


1 


-0.3019395(5)±0.3787735(5)i 


0.87 


-0.1690(2) 


1.1692(2) 


1.1693(3) 


Spin-| 


3 


0.63(l)±0.45(l)i 


1.05 


-1.8(5) 


2.7(5) 


2.4(6) 


1 


0.09477(2)±0.64117(5)i 


0.88 


-0.174(5) 


1.185(5) 


1.185(1) 


2 


-0.0654(5)±0.7113(4)i 


0.97 


-0.18(3) 


1.21(2) 


1.22(3) 


1 


-0.52924(2)±0.33797(2)i 


0.85 


-0.177(5) 


1.184(5) 


1.188(5) 


Spin-2 


2 


-0.842(5) 


1.02 


0.130(4) 


1.2(5) 


0.3(4) 


1 


0.3767(2)±0.6401(l)i 


0.90 


-0.16(3) 


1.19(1) 


1.19(3) 


2 


0.302(6)±G.727(8) 


0.95 




1.3(4) 


1.2(3) 


4 


0.215(15)±0.805(15)i 


1.00 








1 


-0.22561(2)±0.68247(4)i 


0.87 


-0.16(2) 


1.194(6) 


1.192(6) 


2 


-0.394(5)±0.700(6)i 


0.97 




1.8(6) 


1.6(4) 


1 


-0.64890(4)±0.28696(4) 


0.86 


-0.180(5) 


1.197(6) 


1.194(6) 


3 


-0.685(15)±0.485(15)i 


1.01 




2.3(5) 


1.4(3) 


Spin-| 


1 


0.5501(3)±0.5842(2)i 


0.91 


-0.4(1) 


1.19(2) 


1.19(4) 


3 


0.522(5)±0.645(10)i 


0.94 


-1.2(4) 






1 


0.0612(2)±0.7759(2)i 


0.88 


-0.2(1) 


1.20(3) 


1.19(2) 


3 


-0.03(l)±0.83(l)i 


0.94 








1 


-0.4227(l)±0.6400(l)i 


0.87 


-0.20(5) 


1.185(15) 


1.21(3) 


3 


-0.575(5)±0.61(2)i 


0.95 








3 


-0.665(15)±0.53(l)i 


0.97 








1 


-0.7213(2)±0.24595(15)i 


0.87 


-0.175(25) 


1.20(2) 


1.20(2) 


4 


-0.745(15)±0.39(2)i 


0.96 








Spin-3 




-0.92(1) 


1.01 








1 


0.6608(4)±0.5232(5)i 


0.93 




1.20(3) 


1.20(3) 


3 


0.645(15)±0.595(15)i 


0.96 


-1.4(5) 


2.0(5) 


2.0(5) 


1 


0.2729(3)±0.7730(4)i 


0.90 




1.20(4) 


1.19(4) 


4 


0.220(15)±0.840(15)i 


0.95 




1.6(4) 


1.6(4) 


1 


-0.1686(l)±0.7902(l)i 


0.89 


-0.19(3) 


1.20(2) 


1.20(2) 


2 


-0.275(5)±0.825(5)i 


0.95 




1.2(3) 


1.2(3) 


1 


-0.54955(5)±0.58351(3)i 


0.88 


-0.20(4) 


1.196(6) 


1.197(5) 


2 


-0.68(l)±0.54(l)i 


0.95 




1.1(4) 


1.0(4) 


1 


-0.76925(10)±0.21430(5)i 


0.88 


-0.185(25) 


1.205(15) 


1.205(15) 
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The distribution of singularities is shown in figure [1|. A remarkable feature of 
the singularity distribution is its regularity. As S increase the complex singularities 
move closer to the perimeter of the physical disc and the distance between the various 
singularities become more uniform. In the limit S* ^ oo it thus seems likely that the 
singularities will converge onto the unit circle. 

Figure 1. The distribution of singularities in the complex u-plane. In all cases the 
circle has radius Uc- 




We find the very old conjecture by Fox and Guttmann that the number of 
singularities inside the physical disc equals qS — 2, where q is the coordination number 
of the lattice (g = 4 for the square lattice), to be invalid for S* > 1. Recently, Matveev 
and Shrock studied the distribution of zeros of the partition function of the square 
lattice Ising model for S = 1, |, and 2. They conjectured that all divergences of the 
magnetisation occur at endpoints of arcs of zeros protruding into the ferromagnetic 
phase and that there are 4[S'^] — 2 such arcs for S > 1, where [x] denotes the integer 
part of X. Our analysis seems to confirm these conjectures for the magnetisation series 
up to 5* = 2. In particular we find evidence of singularities close to the endpoints located 
by Matveev and Shrock for these spin values. 

The estimate for 7' at the singularity m_ = — 1 of the spin-| susceptibility and the 
estimates for the spin-1 series are based on the low-temperature series we published 
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elsewhere ^ |^. The estimate for 7' of the spin-i case is consistent with the exact value 
7' = I also reported by Matveev and Shrock ITS . 



From table |^ we observe that the exponents at the singularities in the complex 
plane which are well converged (those marked 1) appear to be independent of S. In 
the case of integer spin it appears that the exponents associated with the singularity 
on the negative u-axis equal those at Uc- While the exponents are independent of S, 
note that they do depend on the lattice structure [ITH], so a much weaker version of 



universality holds at the non-physical singularities. In all these cases we observe that 
the Rushbrooke inequality 



a' + 2/3 + 7' > 2, (12) 

is satisfied, and it does indeed seem quite possible that the exponents satisfy the equality 
in Eq. (|12|). At the remaining singularities the errors on the exponent estimates are too 
large to make any such assertion. 



E-mail or www retrieval of series 



The low-temperature series for the spin-S* Ising model can be obtained via e- 
mail by sending a request to iwan@maths.mu.oz.au or via the world wide web on 
http: / / www.maths.mu.oz.au/ ^iwan/| by following the instructions. 
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